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ABSTRACT: The edge crack torsion (ECT) test is designed to initiate mode III 
delamination growth in composite laminates. An ECT specimen is a rectangular 
laminate, containing an edge delamination at the laminate mid-plane. Torsion load 
is applied to the specimens, resulting in relative transverse shear sliding of the 
delaminated faces. The test data reduction schemes are intended to yield initiation 
values of critical mode III strain energy release rate, Guic, that are constant with 
delamination length. The test has undergone several design changes during its 
development. The objective of this paper was to determine the suitability of the 
current ECT test design as a mode III fracture test. To this end, ECT tests were 
conducted on specimens manufactured from IM7/8552 and specimens made from 
S2/8552 tape laminates. Several specimens, each with different delamination 
lengths are tested. Detailed, three-dimensional finite element analyses of the 
specimens were performed. The analysis results were used to calculate the 
distribution of mode I, mode II, and mode III strain energy release rate along the 
delamination front. The results indicated that mode Ill-dominated delamination 
growth would be initiated from the specimen center. However, in specimens of 
both material types, the measured values of Guic exhibited significant dependence 
on delamination length. Furthermore, there was a large amount of scatter in the 
data. Load-displacement response of the specimens exhibited significant deviation 
from linearity before specimen failure. X-radiographs of a sample of specimens 
revealed that damage was initiated in the specimens prior to failure. Further 
inspection of the failure surfaces is required to identify the damage and determine 
that mode III delamination is initiated in the specimens. 



Glossary of Terms 


a Insert length. 

A Constant from relation between stiffness and normalized insert length. 
b ECT specimen width. 

B Delamination front element length in x-direction. 

Cf r Test compliance. 

C spc ECT specimen compliance. 

C su b Sublaminate compliance. 

C sub Average sublaminate compliance. 

Cgys System compliance. 

5 Cross-head displacement. 

Aa Delamination front element length in y-direction. 

AA Virtual area. 

En In-plane modulus along fiber direction. 

E22 In plane modulus perpendicular to fiber direction. 

E33 Transverse modulus. 

G12 In-plane shear modulus. 

G13 Transverse shear modulus. 

G23 Transverse shear modulus 

Gi Mode I strain energy release rate. 

Gn Mode II strain energy release rate. 

Gm Mode III strain energy release rate. 

Gic Critical mode I strain energy release rate (initiation value). 

Guc Critical mode II strain energy release rate (initiation value). 

G/ff C Critical mode III strain energy release rate (initiation value). 

Gt Total strain energy release rate. 

G T Average total strain energy release rate. 
h,4 Thickness of Sublaminate. A. 

hfi Thickness of Sublaminate. B. 

I Distance separating load and support points along ECT specimen length. 
L ECT specimen length. 

m Constant from relation between stiffness and normalized insert length. 

P Applied load. 

pma x M ax j mum test load. 

P C NL Load at onset of non-linearity. 

P c 5% Load at intersection of 5 % offset curve and load displacement response. 
uu Delamination front node displacement in X-direction, column L, row i. 

vu Delamination front node displacement in Y-direction, column L, row i. 

W Distance separating load and support pins along specimen width. 
wu Delamination front node displacement in Z-direction, column L, row i. 
Xu X-direction delamination front force in column L, row i. 

Yu Y-direction delamination front force in column L, row i. 

Zu Z-direction delamination front force in column L, row i. 



Introduction 


Fiber reinforced composite materials are susceptible to delamination due to 
out-of-plane and shear loads, yielding a reduction in structural integrity and 
therefore reduced service lifetime [1]. Delamination in a structural laminate may 
consist of the mode I (opening), mode II (sliding shear) and mode III (scissoring 
shear) components of strain energy release rate. Considerable attention has been 
paid to mode I and mixed mode I/II fracture, resulting in standard testing techniques 
for characterizing such fracture behavior [2-3 respectively]. Mode II delamination 
has received some attention with the development of the four point bend end-notch 
flexure test (4ENF) [4]. Recent work [5], however, indicates there are a number of 
problems with the test method, that need to be resolved prior to standardization. 

Mode III delamination has received considerable attention in the research 
community. Several test techniques have been proposed for characterizing mode III 
delamination fracture, including a split cantilever beam [6] and a crack rail shear 
specimen [7]. A more comprehensive description of developed mode III test 
methods is given in work detailing an anti-clastic plate bend test for mode III 
fracture [8], The most commonly investigated mode III fracture test method is the 
edge crack torsion (ECT) test [9], which is the main subject of this paper. An ECT 
specimen is a rectangular laminate of tape composite material, containing an edge 
delamination at the mid-plane of the specimen. Equal and opposite moment arms 
are applied to the specimen ends (indicated by arrows between w in Fig. 1). The 
applied moment arms generate a torsion load in the specimen, resulting in relative 
shear sliding of the delaminated faces as depicted in Fig. 1. The original specimen 
stacking sequence was [90/(+45/-45) n /(-45/+45) n /90] s , with n=3 or 4 depending on 
the composite material used to manufacture the specimens. The orientation 
corresponds to the coordinate system depicted in Fig. 1. .It is assumed that 
delamination growth takes place along the 90° direction at the 90°/90° ply interface, 
constituting mode III delamination. 

The original ECT test fixture consisted of a load frame whereby the 
specimen was positioned between three support pins, and load was applied via a 
fourth pin, as illustrated in Fig. 2a. A number of studies concluded that this 
specimen produced mode Ill-dominated loading at the delamination front in the 
center of the specimen, away from the loading points [10-11], A round robin 
exercise was organized by the ASTM D30 committee to characterize the test 
method [12]. The ECT tests were undertaken in five independent laboratories on 
specimens manufactured from the toughened carbon/epoxy system, HTA/6376. 
Large scatter in calculated fracture toughness values were reported, and some 
laboratories reported significant nonlinearity in load-displacement response of the 
specimens. Consequently, the test frame was modified, resulting in a symmetrical 
load application via two pins, as illustrated in Fig. 2b. A loading beam was used to 
allow uniform load transfer from test machine to the loading pins. Zero-degree 
plies were added to the specimen, yielding the stacking sequence [90/0/(+45/-45) n /(- 
45/+45)„/0/90] s . 

A second round-robin study involved ECT tests (using the modified 
specimen and test fixture) conducted on specimens manufactured from S2/8552 
glass/epoxy tape and IM7/8552 carbon/epoxy tape [13]. Delaminations were found 



to grow along the 90/0 interface in some glass/epoxy specimens, instead of along 
the intended 90/90 interface [13]. Additionally, it was shown that a significant 
mode II component of strain energy release rate, Gn, was generated near to the 
location of the loading pins [13]. It was argued that increasing the specimen length, 
L, would reduce this mode II component and therefore promote mode III 
delamination in the center portion of the specimen [13]. Consequently, the 
specimen length was increased from 82.5-mm to 108-mm. 

The objective of this paper was to characterize the current ECT test method 
to detennine its suitability for inducing mode III delamination growth. Tests were 
conducted on specimens manufactured from IM7/8552 and specimens made from 
S2/8552 tape laminates. Specimens with insert lengths (normalized by specimen 
width, b), a/b, of 0.2, 0.3, 0.4, 0.5, and 0.6 were tested. Two data reduction 
techniques were used to calculate the critical mode III strain energy release rates, 

G/i/c. 

A selection of specimens was also inspected using a dye-penetrant X-ray 
technique. Detailed, 3D finite element analyses of all specimens were conducted. 
Analysis results were used to calculate distribution of strain energy release rate 
components along the edge delamination front. Findings from the ECT tests and 
analysis were used to determine the suitability of the current data reduction methods 
for calculating Gm c . 


Experimental Procedures 


Specimen and Materials 

An ECT specimen is a rectangular laminate of tape composite material, 
Fig. 1. Dimensions of the specimen are also provided in the figure. A 13pm-thick 
PTFE (Teflon®) insert was positioned at the mid-plane of the specimen to introduce 
an edge delamination crack. Stacking sequence of the specimen was a function of 
material used. Specimens manufactured from IM7/8552 carbon/epoxy tape had the 
stacking sequence, [90/0(+45/-45) 2 /(-45/+45)2/0/90] s . Additional cross-plies were 
added to specimens manufactured from S2/8552 glass/epoxy tape laminate, yielding 
the stacking sequence, [90/0(+45/-45)3/(-45/+45)3/0/90] s . Ply orientations 
correspond to the coordinate system given in Fig. 1. The fiber volume fraction of 
the specimens was 60%. Specimens were cured in an autoclave using cure cycle 
suggested by the composite material manufacturer. All specimens were 
manufactured was by Bell Helicopter Textron. 

Specimens with five different insert lengths were manufactured and tested. 
The five nonnalized insert lengths, a/b, were 0.2, 0.3, 0.4, 0.5 and 0.6. Three repeat 
specimens of each insert length were tested, resulting in a total of 15 IM7/8552 
ECT specimens and 15 S2/8552 ECT. 

ECT Test Fixture 

The load frame used for ECT testing is a symmetric two-point test fixture, 
Fig. 2b. Two support points are located diagonally to each other at the corners of 



the test fixture. Two columns are located at the opposite corners to the support 
points. The columns contain vertical holes lined with spherical bearings. 
Specimens are placed onto the support points. Loading pins are placed through the 
holes in the columns, in order make contact with the specimen. The bearings lining 
the column holes reduce sliding friction between loading pin and the column. A 
loading beam is placed onto the loading pins, as illustrated in Fig. 2b. During an 
ECT test, load is applied at the center of the loading beam, in the direction indicated 
in Fig. 2b. The load is equally distributed to the ECT specimen through the load 
pins. Three guide pins are positioned on the fixture to enable precise specimen 
alignment. Two of these pins are visible in Fig. 2b. The third pin is hidden by the 
left-hand load-pin column. 

ECT Tests 

Prior to testing, all specimens were dried in an oven at 104°C for a 7-day 
period. Specimens were placed in a dessicator after the drying period and were 
tested within 1 day after removal from the oven. 

Specimen dimensions were measured to the nearest 0.05mm. Measurements 
of the width (dimension b ) were taken at the mid-point and at a distance 6.35 mm 
from each end in the y-direction (Fig. 1) resulting in three measurements. 
Measurements for the length (dimension L) were taken at the mid-point and at a 
distance 6.35 mm from each end in the x-direction, again resulting in three 
measurements. Measurement of the thickness were taken at the center of the 
specimen and at the intersections of the lines created a distance 6.35mm from each 
edge resulting in five measurements. All specimen dimensions were taken as the 
average of the corresponding measurements. 

The ECT tests were conducted using a servo-hydraulic test machine. 
Specimens were placed into the ECT test fixture, such that contact was made with 
all three guide pins, ensuring precise specimen alignment. After being leveled, 
specimens were loaded under displacement control at a rate of 1.3mm/min until 
failure. Specimens were unloaded at a rate of 5mm/min. Applied load, P, and 
cross-head displacement, d (referred to as displacement in remainder of paper), 
were recorded during each test using data acquisition software on a computer 
connected to the test machine. A load-displacement response typical from tests on 
both material types is given in Fig. 3. In most specimens, failure was indicated by a 
sudden reduction in load (plot 1 in Fig. 3). The maximum load applied to each 
specimen, P" iax was recorded. Additionally, the load corresponding to the onset of 
nonlinearity of the load-displacement response, P C NL , was calculated using the 
technique detailed in Appendix A. An illustration of the location of P™ ax and P C NL 
on the load-displacement response is given in Fig. 3. IM7/8552 specimens with the 
largest insert length ( a/b = 0.6) exhibited stable failure, as illustrated by plot 2 in 
Fig. 3. In this case, a line with slope 5% less than the original load-displacement 
response was superimposed onto the plot. The load corresponding to the 
intersection of the two curves, P ^ % , was recorded instead of P c '" ax . 



The total specimen compliance, C/ r , following each ECT test was calculated 
by taking the reciprocal of the slope of the load-displacement plot, as illustrated in 
Fig. 3. 

Specimens containing the two smallest insert lengths ( a/b ratios 0.2 and 0.3) 
were split about the specimen mid-plane, creating two sublaminates, labeled A and 
B. The thickness of each sublaminate, Ha and hs, was measured to the nearest 
0.05mm. Each sublaminate was returned to the ECT test fixture and loaded under 
displacement control at a rate of 1 .3mm / min to a load equal to the maximum value 
attained in the original test. Load and displacement were again recorded throughout 
the test. The total compliance of each sublaminate, C su b, was then calculated 
following the same technique used for the intact ECT specimens. An average of the 
compliance of sublaminates A and B, C sub , was then calculated. 

The system compliance (compliance of test machine and ECT load frame), 
C sys , was measured by loading a steel block positioned in the ECT test frame, up to 
the maximum load observed during the ECT testing. Load and displacement were 
recorded and C sys was taken as the reciprocal of the slope of the load-displacement 
response. 

After testing, all ECT specimens were split about the mid-plane and the 
insert length was measured at three locations along the delamination front length. 
The insert length, a, was then taken as the average of the three measurements. 

Two additional carbon/epoxy specimens with normalized insert lengths of 
0.3 were tested. The specimens were loaded to a level higher than the values of 
P C NL observed in previous tests on duplicate specimens, but lower than the 
maximum load, P c max . The specimens were unloaded and held at a constant 
displacement when the load reached approximately half the maximum test load. To 
prepare fore X-ray inspection, a zink-iodide-based dye penetrant was then applied 
to the delamination edges of both specimens, taking care not to spill any penetrant 
on the specimen surfaces. This process lasted approximately 1 minute. The 
specimens were then unloaded. The load-displacement response of the specimens 
was recorded and the technique described in Appendix A was used to confirm that 
P^ L had been reached. 

Data Reduction Methods 

Two data reduction schemes were used to calculate initiation values of the 
critical strain energy release rate, Gm c - The first data reduction method employed a 
multi-specimen compliance calibration procedure. The second technique utilized a 
closed-form solution for mode III strain energy release rate, which was derived 
from laminated plate theory (LPT) in a previous study [10]. 


Compliance Calibration Method 

The compliance of each ECT specimen, C spc , was calculated by subtracting 
the system compliance from the test compliance ( C spc = Cf r -C sys ). The stiffness 
( 1 / C spc ) of all fifteen specimens of each material type were then plotted on the same 



graph as a function of normalized insert length, a/b. Linear regression analysis was 
performed to determine the constants, A and m of the following expression for 
specimen stiffness [9]: 


-!—= /[l-rr(a/bj\ (1) 

spc 

Only the constant, m, was used in the data reduction. 

The perceived critical strain energy release rate of each EOT specimen was 
calculated based on the maximum critical load, P™ ax , and the load corresponding to 
the onset of nonlinearity, P C NL using the following expressions [9]: 


Gfiic ( max ± 


mcjpr[ 

2! 6(1 - rr(a/b )] 


( 2 ) 


G£(NL)- 


2 1 6(1 - rr(a bj\ 


( 3 ) 


The superscript on the left hand side of Eqns. 2 and 3 denote compliance 
calibration. The parameter, /, is the distance separating the load pins along the 
specimen length and b is the specimen width, Fig. 1. 

The perceived critical strain energy release rates, G^ c (max and G(( c ( NL ), 
were then plotted as functions of normalized insert length, a/b. This was repeated 
for both materials tested. It is noted that the resulting critical strain energy release 
rate values correspond only to fracture initiation. Hence, fracture resistance effects 
such as fiber bridging should not be present. The values of G^ c (max and 
G((( c (NL) should therefore be independent of insert length. 


Laminated Plate Theory (LPT) Method 


The perceived critical strain energy release rate of specimens with 
normalized insert lengths, a/b, of 0.2 and 0.3 were calculated using the LPT 
method. As with the compliance calibration technique, perceived critical strain 
energy release rates corresponding to P c m3X and P C NL were calculated for each 
specimen using the following expressions [10]: 
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The parameter L is the total ECT specimen length and W is the distance separating 
the load pins along the specimen width, Fig. 1. 

The values of G(( T (max and G((J ( NL ) were then plotted as functions of 
normalized insert length, a/b. This was repeated for both materials tested. The 
values were superimposed onto the plots of G^ c (max and G^ c (NL) versus a/b, 
providing a comparison of the values calculated using the two data reduction 
methods. 

Dye-Penetrant X-Ray Imaging 

X-radiographs were taken of the two carbon/epoxy ECT specimens, 
penetrated with ink-based dye, using a Pantak Seifert X-ray system. This was done 
to determine whether damage was initiated from the insert after P C NL but before the 
maximum load, P™ ax , was reached. Radiographs were taken along three sections of 
each specimen. The images were then stitched together using graphics post 
processing software, yielding complete images of both specimens. 


Numerical Analysis 


Finite Element Models 

Three-dimensional finite element models were constructed of the IM7/8552 
and S2/8552 ECT specimens. Specimens with each nonnalized insert length were 
modeled. A summary of the specimen dimensions is given in Fig. 1. The models 
were constructed using the commercial code, ABAQUS version 6.3. Solid, eight- 
node brick elements were used to represent the specimens. A composite layer 
option was used to represent specimen stacking sequence, whereby one layer of 
elements was used to represent one or more plies. In this case, the orthotropic ply 
properties were oriented according to the specimen stacking sequence. An image of 
a finite element model (showing displaced geometry) is given in Fig. 4, illustrating 
the stacking sequence of a carbon/epoxy specimen. The edge delamination was 
modeled by including elements with coincident nodes on the plane of the 
delamination. A fine mesh was used in the vicinity of the delamination front to 
accommodate for the rapid change in strain field. The element thickness at the 
delamination front (in the y and z-axes) was one ply thickness as illustrated in Fig. 
4. A similar meshing technique was adopted during an analysis of a double 
cantilever beam specimen [14]. Contact elements were used at the edge 
delamination plane to prevent mesh interpenetration during the analysis. Relative 
sliding between points in the delamination region was assumed frictionless. A 
prescribed displacement in the z-axis was applied at nodes corresponding to the 
point of contact of the loading pins. The same displacement was prescribed to each 
model, therefore simulating displacement control used during actual ECT tests. 
Displacement values were chosen to ensure an elastic response from the specimens, 
which was 2 -mm in all cases. Nodes positioned at the locations of the support pins 
were constrained from displacement in the z-axis to represent contact between 



specimen and the pins. Nodes positioned at the locations where contact takes place 
between the specimen and the load frame guide pins, were constrained from 
displacement in the corresponding axes, as illustrated in Fig. 4. These boundary 
conditions acted to prevent rigid body motion during an analysis run. The 
orthotropic ply properties used to represent IM7/8552 and S2/8552 are presented in 
Table 1. Geometrically nonlinear analyses were performed to facilitate contact in 
the models. 

The reaction loads at the nodes to which displacement was prescribed was 
calculated, and specimen compliance was then determined by dividing prescribed 
displacement by the sum of these reaction loads. 

Virtual Crack Closure Technique 


The virtual-crack-closure-technique (VCCT) [16] was used to calculate the 
mode I, mode II and mode III components of strain energy release rate (G/, Gu and 
Gin respectively) along the delamination front in each finite element model. A 
previously developed ABAQUS user subroutine was used to perform the VCCT 
calculations [17]. The technique works on the principle that the change in stored 
elastic strain energy associated with a small extension of crack area is equal to the 
work done required to close the crack to its original length. In terms of a finite 
element model constructed from solid, 8-node brick elements, G/, Gu and Gui were 
calculated using the following equations [18]: 
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Figure 5 contains an illustration of the delamination front elements typical in the 
finite element models of the ECT specimens. The area AA = Aa B as shown in 
Fig. 5, where AA is the virtual area closed, Aa is the length of the elements at the 
delamination front and B is their width. The subscript in Eqns. 6-8 denote rows 
and columns of nodes as seen in the top view of the delamination front elements in 
Fig. 5b. Capital letters indicate columns and lower case letters indicate rows. 
Hence, Xu, Yu and Z/_, denote the forces at the delamination front in row i, column 
L. The corresponding displacements behind the delamination front at the top face 
of node row X, column L in the x, y and z axes are denoted by un t , v Lx and wlx 
respectively. The displacements at the bottom face of node row X, column L are 
denoted by ul x *, vl-,* and All the forces and displacements are obtained from 
the finite element analyses with respect to the global coordinate system (x, y, z). As 
each analysis was geometrically nonlinear, the forces and displacements were 
resolved into the local coordinate system (x’,y’,z’) using the technique detailed in 



[17]. Equations 6-8 were then used to calculate Gi, Gu and Gm at every node 
located along the delamination front. The strain energy release rate values were 
then plotted as a function of location along the delamination front, x/L. 

The total strain energy release rate at any location along the delamination 
front, Gt, was calculated as the sum of the individual strain energy release rate 
components: 

G t =G, + Gu + Gm (9) 

The total average strain energy release rate across the entire delamiantion length, 
G t , was computed as the integral of the total strain energy release rate divided by 
the delamination length. 


Results and Discussion 


Numerical Analysis Results 

Finite Element Model Verification 

The compliance of each ECT specimen, C spc , was estimated from the finite 
element models and plotted as a function of normalized insert length. Figure 6 
presents specimen stiffness versus normalized delamination length, a/b, calculated 
from the finite element models of the IM7/8552 ECT specimens. Included in the 
figure is a plot of specimen stiffness measured from the corresponding ECT tests. 
The constants A and m of Eqn. 1 were calculated using the finite element analyses 
and experimental data. Analysis and measured stiffness values agree to within 5%, 
indicating the analyses accurately captured the elastic response of the specimens. A 
similar comparison was made from the finite element analyses of the S2/8552 ECT 
specimens. Stiffness versus normalized insert length of these specimens is also 
given in Fig. 6. Again, plots from analysis and experiment are included for 
comparison. For the shorter insert lengths, a/b = 0.2 and 0.3, the agreement 
between analysis and experiment is only within 10%, however, the agreement 
improves dramatically for the largest three insert lengths. Overall, the finite 
element models accurately captured the elastic response of the S2/8552 specimens. 

Strain Energy Release Rate Distribution 

Figure 7a contains plots of Gn and Gm versus distance along the 
delamination front computed from analyses of the IM7/8552 ECT specimens. The 
plots correspond to specimens with the smallest and largest normalized insert 
lengths, a/b = 0.2 and a/b =0.6 respectively. The mode I strain energy release rate 
was found to be negligible in comparison to Gu and Gm for all specimens modeled, 
and is therefore not included in the plots of Fig. 7. For a given location along the 
delamination front, Gu and Gm, were found to decrease with increasing insert 
length. The strain energy release rate distribution is very similar to that reported 
from an analysis of the original ECT specimen [19] (Fig. 2a), with the exception 
that in the current analyses, the distribution was found to be symmetrical about the 



specimen mid-length. The parameter, Gn, peaks at the locations of the load and 
support pins and Gui peaks along the center of the specimen. The load and support 
pins produce a moment arm that cause relative sliding of the delaminated sections 
of the specimen, parallel to the intended direction of delamination growth. 
Consequently, Gn is expected to peak at the location of the load pins. The finite 
element analyses were also used to calculate contact pressure in the region of the 
delamination. Contact pressure was negligible across most of the region, except at 
positions local to the load and support pins. Relative opening of the contact faces 
was also negligible. This supports the finding that Gi is negligible. It should be 
noted, however, that sliding friction of the delaminated surfaces was not considered. 
It was assumed that friction would not influence the strain energy release rate 
distribution along the delamination front. 

In the case of all the IM7/8552 specimens, the peak value of Gm was 
approximately eight times the peak Gn value. Assuming Gni c ^ Guc, the above 
findings indicate that mode Ill-dominated delamination will initiate from the center 
of the specimen insert front. 

Similar trends were found from analyses of S2/8552 ECT specimens, with 
the exception that Gn exhibited a slight increase with insert length, as illustrated in 
Fig. 7b. Again strain energy release rate distribution from analyses of specimens 
with smallest and largest insert lengths are presented in Fig. 7b. The results 
indicate that mode Ill-dominated delamination growth should take place from the 
center region of the specimen. 

Effect of Insert Length on Average Total Strain Energy Release Rate 

Figure 8 contains plots of average total strain energy release rate, G r , 
versus normalized insert length, calculated from analyses of the IM7/8552 and 
S2/8552 ECT specimens. In the case of the IM7/8552 specimens, G T decreases 
with insert length. Considering that the same displacement was prescribed in each 
analysis, the reduction in G T with insert length implies that delamination growth 
will be stable if testing is performed under displacement control. This result 
suggests that the use of Pf ax for calculating Gnic may be inaccurate, as 
delamination may occur prior to the sudden load drop that was observed during 
ECT tests. A similar trend was predicted for the S2/8552 specimens, except that an 
initial increase in G T was predicted for an increase in nonnalized insert length from 
0.2 to 0.3, as shown in Fig. 8. 

Experimental Results 

ECT Specimen Load-Displacement Response 

A sample of experimentally-measured load-displacement responses from 
tests on IM7/8552 and S2/8552 specimens is presented in Figs. 9a and 9b 
respectively. The responses from specimens of the two material type were similar. 
A typical load-displacement plot consisted of an initial nonlinear response, 
corresponding to seating of the load and support pins on the specimens. This was 



proceeded by a linear response, after which the load-displacement response 
deviated from linearity. In most cases, the load reached a peak value, followed by a 
sudden load drop that was assumed to correspond to specimen failure. The 
IM7/8552 specimens containing the largest insert length ( a/b = 0.6) exhibited stable 
failure, shown by the gradual change in slope at the upper portion of the load- 
displacement response. In the case of the IM7/8552 specimens, the deviation from 
linearity, P C NL , ranged from 90% to 93% of the maximum load, Pf ax . The extent 
of deviation from linearity was found to increase with increasing insert length. The 
glass/epoxy specimens exhibited more exaggerated deviation from linearity, with 
P C NL ranging from 50% to 75% of Pf ax . Again the extent of deviation from 
linearity increased with insert length. 

Critical Mode III Strain Energy Release Rate of IM7/8552 Specimens 

The perceived critical mode III strain energy release rate was calculated 
using the compliance calibration method. A plot of specimen stiffness versus 
normalized insert length from tests on both material types is presented in Fig. 10. 
The plots of l/C S p C versus a/b from specimens of both materials exhibit a linear fit, 
validating the relationship of Eqn. 1 . The values of perceived critical mode III 
strain energy release rate were calculated using Eqns. 2 and 3. Figure 11a contains 
the plots of Gjj£ (max and G^(NL) versus normalized insert length for the 
IM7/8552 specimens. The values of G^ c (max and G'jj^NL) are in good 
agreement for the two smallest insert lengths. However, G ^ (max begins to 
deviate from Gf„ c (NL ) when a/b > 0.3, where the values of G^ (max begin to 
increase with a/b while G^^NL) remains relatively constant. The G^j^NL) 
values exhibited significant scatter. Similar scatter was observed in data generated 
from ECT tests conducted on glass/epoxy specimens [20]. The findings imply that 
the values of G ^ (max do not correspond to Gnic (of IM7/8552) over the largest 
three insert lengths. In theory, the values are representative of fracture initiation, 
and therefore the measured critical mode III strain energy release rate should be 
constant with insert length. 

The perceived critical mode III strain energy release rates calculated using 
the laminated plate theory method (Eqns. 3 and 4) are also plotted in Fig. 11a. The 
values of G[j( c T (NL ) agree very well with the G ^ (NL ) values for specimens with 
the two smallest insert lengths. Similarly, values of G^J (max are in good 
agreement with Gf£(max. 

Averages of G^ (NL ) and G^ (max were calculated and plotted on the 
column chart presented in Fig. lib. An average value for the critical mode II strain 
energy release rate of IM7/8552 is also given in the plot for comparison. The mode 
II fracture data was measured from four point bend end-notch flexure tests 
conducted on IM7/8552 specimens [21]. Error bars correspond to one standard 
deviation above and below the average values represented by the columns. The 
average values are also printed at the top of each column. The column representing 



Guc is an average of only three data points, so caution is taken as to the statistical 
significance of the value. In all cases, scatter in the data is high, resulting in overlap 
of the (max and Gn c values. The average value of G^ (NL ) is 17% lower 
than the average Gn c . It is expected that Gnic should be equal to or greater than Gn c . 
Hence, even though G^ (NL ) appears consistent with insert length, this value also 
may not be representative of Gm c . 


Critical Mode III Strain Energy Release Rate of S2/8552 Specimens 


Figure 12a contains plots of G^fNL), Gf c (max, 


'CC 


and 


f' LPT , 

b,, 


. m c (max measured from tests on the S2/8552 ECT specimens. Again, the values 
of G^ (max increase with insert length, indicating the data is not representative of 
Gmc . The values of Gf^ c (NL ) remain relatively constant with insert length. There 
is a greater difference between G^(NL) and G^ (max in comparison to the 


IM7/8552 values. This is because the glass/epoxy specimens exhibited significantly 
more deviation from linearity before failure. There is also a significant difference 
between fracture data calculated using the laminated plate theory method and those 
values found from compliance calibration. Reasons for this difference are not 
understood, although they may be associated with the nonlinear specimen response 
during loading. 

Average values of G^ (NL ) and G^ (max are presented in the column 


chart in Fig. 12b. Values of Gn c for S2/8552 [21] are also contained in the chart. 
Again, the mode II value is an average of only three data points. The results show 
that G(jK (max is much greater than Gn c . Furthermore, the average of G^' (NL ) is 
20% lower than Gu c , suggesting that the G^fNL) values may not be 
representative of Gm c . 


X-Radiograph Imaging 


Figure 13a contains X-radiographs taken of the two carbon/epoxy ECT 
specimens (Specimens A and B) that were loaded to a level between P C NL and Pf ax . 
The load-displacement response of both specimens is shown below the 
X-radiographs in Fig. 13b. Regression analysis was used to extrapolate the linear 
portion of the load-displacement responses, and the method detailed in Appendix A 
was used to determine the values of P C NL . It is seen that both specimens exhibited 
some deviation from linearity. The corresponding X-radiographs show that some 
type of damage initiated from the insert fronts (indicated by the light patches on the 
surfaces). The damage is likely to consist of delamination and splitting of the 90° 
plies adjacent to the insert plane. The findings suggest that P C NL corresponds to the 
initiation of these energy-absorbing mechanisms. Therefore, the values of 
Gf£ (max and G(uJ (max do not likely correspond to delamination initiation. 
Additionally, the increase in extent of deviation from linearity with insert length 
explains the observed dependence of G^ (max with a/b. 


It was observed that the average values of Gf£ (NL ) for each material was 
less than the corresponding mode II values. This indicates that the damage 
observed in the X-radiographs may not constitute initiation of mode III 
delamination growth. Detailed fractographic analysis is required to confirm this 
hypothesis. 


Discussion 

The numerical and experimental results appear to indicate that the fracture data 
generated from ECT tests on the IM7/8552 and S2/8552 specimens do not 
correspond to mode III fracture initiation. All specimens exhibited some deviation 
from linearity, indicating that damage growth may be stable. This was corroborated 
by the findings from the finite element analysis, showing that the average total 
strain energy release rate decreases with insert length, Fig. 8. 

Data from finite element analyses of the ECT specimens also indicated, that 
mode Ill-dominated delamination growth should initiate from the center of the 
specimens. 

X-radiograph images taken of specimens loaded between P C NL and P c max 
indicated that some damage takes place before the maximum test load, P™ ax , is 
reached. Detailed fractographic analyses of the specimen failure surfaces is needed 
to identify the type of damage revealed by the X-radiographs. 

In summary, from the results of studies conducted in the present work, the 
current data reduction schemes may not be suitable for calculating critical mode III 
strain energy release rate. Furthermore, detailed inspection of the specimen failure 
surfaces is required to determine whether mode III delamination was initiated in the 
specimens. 


Concluding Remarks 

Edge crack torsion tests were conducted on specimens manufactured from 
IM7/8552 and S2/8552 tape laminates. In most instances, failure was indicated by a 
sudden drop in load on the load-displacement response. Two data reduction 
schemes (compliance calibration and laminated plate theory) were used to calculate 
perceived critical mode III strain energy release rates. Values calculated using the 
maximum test load ( G(j£ (max) exhibited significant dependence on specimen 
insert length. As the values are meant to represent fracture initiation, the 
dependence of G(jK(max with a/b should not be present. It was therefore 
concluded that the values do not represent initiation of mode III delamination 
growth. The Gf£ (NL ) values were moderately insensitive to a/b, however, the 
average value was less than the average mode II value measured in previous work 
[21]. This implied that these values also may not reflect initiation of mode III 
delamination growth. 



Two IM7/8552 specimens both with normalized insert lengths of 0.3 were 
tested. The specimens were loaded to a level above P C NL but were not allowed to 
fail (load did not reach P™ ax ). X-radiograph images of the specimens revealed that 
damage took place. The observed deviation from linearity of the specimen load- 
displacement responses suggest the damage growth was stable. This was supported 
by the findings from the finite element analyses, that average total strain energy 
release rate decreases with insert length for a constant specimen displacement. The 
analyses also indicated that specimens should exhibit mode III dominated 
delamination growth about the specimen center. Detailed inspection of the 
specimen failure surfaces is required to confirm the numerical analysis findings. 
This is also necessary to identify the appropriate data reduction strategies if 
delamination growth is found to be mode III dominated. 
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Appendix A 

Calculating Load at Deviation from Linearity 

The load-displacement response from an ECT test on a glass/epoxy specimen with a 
normalized insert length of 0.2 is shown in Fig. Al. The response is initially 
nonlinear, corresponding to seating of the load and support pins on the ECT 
specimen. This is proceeded by a linear section, corresponding to the elastic 
response of the specimen. At some load, P C NL , the response begins to deviate from 
linearity until a maximum load is reached, followed by a sudden load drop, 
indicating specimen failure. A linear regression analysis is performed to estimate 
the linear fit corresponding to the linear portion of the load-displacement response 
(load range bounded by broken lines in Fig. Al). The load is calculated using the 
linear fit over the entire displacement range of the load-displacement response. A 
plot of this linear fit load, Plr, is superimposed onto the original load-displacement 
plot in Fig. Al. The difference between the linear fit load and the original 
experimental load (Plr-Pexp) is plotted as a function of displacement, as shown in 
Fig. A2. The onset of deviation from linearity was then determined by recording the 
displacement value after which the linear regression and experimental load values 
begin to diverge, as shown in Figure A2. The load at deviation from linearity, P C NL , 
is then found from the original load-displacement plot, as illustrated in Fig. Al. 



TABLE 1. MATERIAL PROPERTIES [15]. 


S2/8552 Unidirectional Glass-Epoxy Prepreg 

E n =47.71 GPa 

E 2 2 = 12.27 GPa 

£33 = 12.27 GPa 

Vi 2 = 0.278 

v i3 = 0.278 

V 23 = 0.403 

G 12 = 4.83 GPa 

Gi 3 = 4.83 GPa 

G 23 = 4.48 GPa 

IM7/8552 Unidirectional Carbon-Epoxy Prepreg 

£"11 = 161.0 GPa 

£ 2 2 = 11-38 GPa 

£33 = 11.38 GPa 

V i2 = 0.32 

V13 = 0.32 

V 23 = 0.436 

Gi 2 = 5.17 GPa 

Gi 3 = 5.17 GPa 

G 23 = 3.98 GPa 


h A +h B 


L=108mm, /=76mm 
b=38mm, W=32mm 
h A +h B =4.5mm (IM7/8552) 
h A +h B =7.5mm (S2/8552) 


0° = Principle fiber direction. 



FIGURE 1. Schematic of ECT specimen. 
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FIGURE 2. (a) Original ECT test fixture, (b) Current ECT test fixture (showing 

deformed specimen). 




FIGURE 3. Example load displacement responses from ECT tests. 



FIGURE 4. Finite element mesh of IM7/8552 ECT specimen with a/b= 0.6. 



FIGURE 5. Delamination front region of finite element models. 



Experiment (IM7/8552) 
Analysis (IM7/8552) 
Experiment (S2/8552) 
Analysis (S2/8552) 


IM7/8552 

A=1 160N/mm, m=1.01 [Exp] 

A=1 131N/mm, m=0.953 [Analysis] 

S2/8552 

A=2658N/mm, m=0.968 [Exp] 
A=2419N/mm, m=0.886 [Analysis] 


FIGURE 6. Specimen stiffness versus a/b. Comparison of experiment and analysis. 



FIGURE 7. Computed strain energy release rate distribution across delamination 
front (a) IM7/8552 specimens, (b) S2/8552 specimens. 
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FIGURE 8. Computed average total strain energy release rate versus a/b. 



FIGURE 9. Load-displacement responses from ECT tests (a) IM7/8552 specimens, 

(b) S2/8552 specimens. 



IM7/8552 

A=1 144N/mm, m=0.984 

S2/8552 

A=2602N/mm, m=0.953 


FIGURE 10. Stiffness versus a/b from ECT tests. 



FIGURE 1 1 (a) Perceived critical mode III strain energy release rate versus insert 
length of IM7/8552. (b) Average perceived critical strain energy release rate values. 



FIGURE 12 (a) Perceived critical mode III strain energy release rate versus insert 
length of S2/8552. (b) Average perceived critical strain energy release rate values. 
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FIGURE 13 (a) X-radiographs of IM7/8552 specimens, (b) Load displacement 

response of the specimens. 
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FIGURE Al. Load-displacement response typical from an ECT test. 




